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The effect of Coulomb interaction in an ensemble of Dirac fermions on the formation of superconducting 
pairing in monolayer and bilayer doped graphene is studied using the Kohn-Luttinger mechanism disregarding 
the Van der Waals potential of the substrate and impurities. The electronic structure of graphene is described 
using the Shubin-Vonsovsky model taking into account the intratomic, interatomic, and interlayer (in the case 
of bilayer graphene) Coulomb interactions between electrons. The Cooper instability is determined by solving 
the Bethe-Saltpeter integral equation. The renormalized scattering amplitude is obtained with allowance for 
the Kohn-Luttinger polarization contributions up to the second order of perturbation theory in the Coulomb 
interaction. It plays the role of effective interaction in the Bethe-Salpeter integral equation. It is shown that 
the allowance for the Kohn-Luttinger renormalizations as well as intersite Coulomb interaction noticeably 
affects the competition between the superconducting phases with the /—wave and d + id—wave symmetries 
of the order parameter. It is demonstrated that the superconducting transition temperature for an idealized 
graphene bilayer with significant interlayer Coulomb interaction between electrons is noticeably higher than 
in the monolayer case. 


1. INTRODUCTION 

Graphene is of considerable interest for fundamental 
physics and for applications due to its peculiar trans¬ 
port, pseudorelativistic, and quantum-electrodynamic 
properties m- This combination of graphene proper¬ 
ties is primarily determined by its unique gapless energy 
structure consisting of cone-shaped valence and conduc¬ 
tion bands contacting at the corners of the first Brillouin 
zone (Dirac points) [4]. It has been established that 
electrons propagating in graphene near Dirac points re¬ 
semble massless fermions with linear dispersion and 
are characterized by the minimal conductivity for a zero 
charge carrier concentration EH, a high mobility [ZH3, 
Klein tunneling [TUlITT]. oscillating motion (Zitterbewe- 
gung) [THIISj, universal absorption of light [13], and 
many other properties having no analogs in other phys¬ 
ical systems. 

When in contact with superconductors, graphene ex¬ 
hibits exotic superconducting properties m- In spite 
of the fact that the evolution of the Cooper instability 
in graphene itself has not yet been confirmed, experi¬ 
mental evidence (MS] that graphene in contact with 
conventional superconductors exhibits superconducting 
properties have been obtained. The fact that short 
graphene samples placed between superconducting con¬ 
tacts can be used to construct Josephson junctions in¬ 
dicates that Cooper pairs can propagate coherently in 
graphene. In this connection, it would be interesting to 


find out whether it is possible to modify graphene struc¬ 
turally or chemically to convert it into a magnet m or 
even into a real superconductor. 

It is known theoretically that a model with conic 
dispersion requires the minimal intensity of the pairing 
interaction for the development of the Cooper instabil¬ 
ity [23j . In this connection, there have been several 
attempts to analyze theoretically possible achievement 
of the superconducting state in doped monolayer, as 
well as bilayer, graphene. The role of topological de¬ 
fects in achieving Cooper pairing in a graphene mono- 
layer was studied in [53] . In [^ , the phase diagram for 
spin-singlet superconductivity in a monolayer was con¬ 
structed by Uchoa and Castro Neto in the mean field 
approximation, and the plasmon mechanism of super¬ 
conductivity leading to low superconducting transition 
temperatures in the s-wave channel was studied for real¬ 
istic electron concentration values. The possibility of in¬ 
ducing superconductivity in a graphene monolayer due 
to electron correlations was investigated in ESHH]. 

The situation in which the Fermi level is near one of 
the van Hove singularities in the density of states of a 
graphene monolayer was considered in [28] . It is well 
known that these singularities can enhance magnetic 
and superconducting fluctuations m- In accordance 
with the scenario described in [55] , the Cooper instabil¬ 
ity occurs due to strong anisotropy of the Fermi contour 
for van Hove filling riyH-, which in fact originates from 
the Kohn-Luttinger mechanism m proposed in 1965 
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and assuming the occurrence of superconducting pair¬ 
ing in systems with purely repulsive interaction |31H33] . 
It was noted in |28) that this mechanism can occur in 
graphene because the electron-electron scattering be¬ 
comes strongly anisotropic; for this reason, a channel 
with attraction can be formed when there is a projec¬ 
tion onto harmonics with a nontrivial angular depen¬ 
dence on the Fermi surface. According to the result 
obtained in m, such the Cooper instability in an ideal¬ 
ized graphene monolayer evolves predominantly in the 
d—wave channel and can be responsible for supercon¬ 
ducting transition temperatures up to Tc ^ 10 K de¬ 
pending on the proximity of the chemical potential level 
to the van Hove singularity. An analogous conclusion 
was drawn in |34] , where the Kohn-Luttinger supercon¬ 
ductivity in the vicinity of the van Hove singularity in 
the graphene monolayer was studied by the renormal¬ 
ization group method. 

The possibility of the competition and coexistence 
of the Pomeranchuk instability and the Kohn-Luttinger 
superconducting instability in a graphene monolayer 
was considered in m- In [ 35 ] , it was found in 
experiments with a strongly doped monolayer using 
angle-resolved photoemission spectroscopy (ARPES) 
that multiparticle interactions substantially deform the 
Fermi surface, leading to an extended van Hove singu¬ 
larity at point M of the hexagonal Brillouin zone. The 
features of the ground state were investigated theoret¬ 
ically, and the competition between the ferromagnetic 
and superconducting instabilities was analyzed. It was 
shown that in this competition, the tendency to su¬ 
perconductivity due to strong modulation of the effec¬ 
tive interaction along the Fermi contour (i.e., due to 
electron-electron interactions alone) prevails. The su¬ 
perconducting instability evolves predominantly in the 
/—wave channel in this case [36] . The competition 
between the superconducting phase and the spin den¬ 
sity wave phase at van Hove filling and near it in the 
monolayer was analyzed in |37j by the functionalization 
renormalization group method. It was found that for 
the band structure parameters and the Coulomb inter¬ 
actions obtained by ab initio calculations for graphene 
and graphite monolayers [53], superconductivity with 
the d + id—wave type of symmetry of the order parame¬ 
ter prevails in a large domain near the van Hove singu¬ 
larity, and a change in the related parameters may lead 
to a transition to the phase of the spin density wave. 
According to far away from the van Hove singu¬ 
larity, the long-range Coulomb interactions change the 
form of the d + id—wave function of a Cooper pair and 
can facilitate superconductivity with the /—wave sym¬ 
metry of the order parameter. 


In accordance with the results obtained in |39| . in 
the case of bilayer graphene, ferromagnetic instability in 
the vicinity of the van Hove singularities dominates over 
the Kohn-Luttinger superconductivity. However, the 
Coulomb interaction screening function in the bilayer 
was calculated earlier in m in the random phase ap¬ 
proximation (RPA) in the doped and undoped regimes. 
It was found that the static polarization operator of the 
doped bilayer contains the Kohn anomaly much larger 
than in the case of a monolayer or a 2D electron gas. It 
is well known that the singular part of the polarization 
operator or the Kohn anomaly [4TH43] facilitates effec¬ 
tive attraction between two particles, ensuring a con¬ 
tribution that always exceeds the repulsive contribution 
associated with the regular part of the polarization op¬ 
erator for the orbital angular momenta Z 0 of the 
pair [ 30 ]. Thus, we can expect that the superconduct¬ 
ing transition temperature Tc in the idealized bilayer 
may exceed the corresponding value for the idealized 
monolayer. 

In addition, it was shown in earlier publications |441 
135 ] that the value of Tc can be increased via the Kohn- 
Luttinger mechanism even for low concentrations of 
charge carriers if we consider the spin-polarized or two- 
band situation, as well as a multilayer system. In this 
case, the role of the pairing spins up is played by elec¬ 
trons of the hrst band (layer), while the role of the 
screening spins down is played by electrons of the sec¬ 
ond band (layer). Coupling between the electrons of the 
two bands occurs via interband (interlayer) Coulomb in¬ 
teraction. As a result, the following superconductivity 
mechanism becomes possible: electrons of one species 
form a Cooper pair by polarizing electrons of another 
species [ 111135 ]. This mechanism of superconductivity 
is also effective in quasi-two-dimensional systems. Note 
that odd-momentum pairing and superconductivity in 
vertical graphene heterostructures made up by graphene 
layers separated by boron nitride spaces was considered 
recently by Guinea and Uchoa [ 25 ] . 

In this work, we investigate the Kohn-Luttinger 
Cooper instability in an idealized monolayer and bi¬ 
layer of doped graphene using the Born weak-coupling 
approximation and taking into account the Coulomb re¬ 
pulsion between electrons of the same and of the nearest 
carbon atoms in a monolayer, as well as the interlayer 
Coulomb repulsion in the case of the bilayer. 

The necessity of including the long-range Coulomb 
interaction in calculating the physical characteristics 
of graphene is dictated by the results obtained in 
[ 35 ], where the partly screened frequency-dependent 
Coulomb interaction was calculated ab initio in con¬ 
structing the effective multiparticle model of graphene 
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and graphite. It was found that the intra-atomic re¬ 
pulsion potential in graphene is U = 9.3 eV (an anal¬ 
ogous estimate is given in |46j 1 . which contradicts the 
intuitively expected small value of U and weak-coupling 
limit U < W. The calculations performed in [38] 
have demonstrated the fundamental importance of tak¬ 
ing into account the nonlocal Coulomb interaction in 
graphene: the Coulomb repulsion of electrons at neigh¬ 
boring sites amounts to 14 = 5.5 eV. It should be 
noted that the values of U and V estimated by other 
researchers (see, e.g., [JT]) are much smaller. 


2. IDEALIZED GRAPHENE MONOLAYER 

In the hexagonal lattice of graphene, each unit cell 
contains two carbon atoms; therefore, the entire lat¬ 
tice can be divided into two sublattices A and B. In 
the Shubin-Vonsovsky (extended Hubbard) model, the 
Hamiltonian for the graphene monolayer taking into ac¬ 
count electron hoppings between the nearest and next- 
to-nearest atoms, as well as the Coulomb repulsion be¬ 
tween electrons of the same atom and of adjacent atoms 
in the Wannier representation, has the form 



fSa- 


Here, operators a'j:^{afa) create (annihilate) an electron 
with spin projection ct = ±1/2 at site / of sublattice A; 

= ajf^a/cr denotes the operators of the number of 
fermions at the / site of sublattice A (analogous nota¬ 
tion is used for sublattice B). Vector 5 connects the 
nearest atoms of the hexagonal lattice. In the Hamil¬ 
tonian, the symbol (( )) indicates that summation is 
carried out only over next-to-nearest neighbors. 

Passing to the momentum space and performing the 
Bogoliubov transformation, 

ccikcr = ■iCii(k)ak<T ± Wj2(k)6k<T, * = 1,2, (4) 


2 

Ho — ^-ikO^iko-'^ikcr- (5) 

i—1 kcr 

The two-band energy spectrum is described by the ex¬ 
pressions m 


Ewa = tl|Mk| — 12 /k, ±^2k = —ll|**k| — l2/k, 

where the following notation has been introduced: 

/k = 2cos{V3ky) ±4cos^^fcy^ cos^fca,^. 


Uk 


= E 


g*ki5 _ ^-ikx 


Vs, 


cos I —ky I, 


l^kl = \/3 ± /k. 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


The coefficients of the Bogoliubov transformation have 
the form 

u;i,i(k) = u;;2(k) = Erk, (10) 

V2 |Mk| 


Wl2(k) = -W2l(k) = - 


V2- 


In the Bogoliubov representation, interaction oper¬ 
ator m is defined by the following expression including 
operators aiko- and a 2 ka- ctika ct 2 ka 


Hint — E/ PI*!’®)‘^Lcr'^jpo-4>:iqCTQ:mscr X 

ijlm 
kpqscr 

X A(k + p - q - s) + (11) 

ijlm 

kpqs 

X A(k ± p - q - s). 


where A is the Kronecker symbol, while r||^.;^(k, p|q, s) 
and r^.;^(k, p|q, s) are the initial amplitudes. The 
quantity 

rl';/™(k,p|q,s) = i(Vij;j™(k,p|q,s) ± 

± Hj*;m/(p,k|s,q)), (12) 

l/y/m(k,p|q,s) = Vuq-pWii{'k)wj 2 ip)wi 2 {q)'w^^{s) 

(13) 


describes the intensity of the interaction of fermions 
with parallel spin projections, while the quantity 

ri|);im(k,p|q,s) = C/y;i^(k,p|q,s) ± 

“1“ (P; k|s, q), (14) 

(k,p|q,s) = I/('u;ii(k)u;ji(p)'u;;i(q)u;;;ji(s) ± 


we diagonalize Hamiltonian Hq , which acquires the form 


+Wi2 0i)Wj2 {p)Wi2 iq)w*„,2 (s)) 


(15) 
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Fig. 1. First- and second-order diagrams for the effective interaction of electrons in graphene monolayer and bilayer. 
Solid lines with light (dark) arrows correspond to the Greens functions for electrons with spin projections -fi (“§) 
and energies corresponding to graphene energy bands aj, ai and a^- In diagrams for the monolayer (Section 2), 
subscripts i = j — 1', subscripts I and m can acquire values of 1 or 2. In the case of the bilayer (Section 3), subscripts i 
and j acquire values of 1 or 2, while subscripts I and m acquire values of 1, 2, 3, or 4. Momenta qi and q 2 are defined 
by relations (O. 


corresponds to the interaction of Fermi quasiparticles 
with antiparallel spin projections. Indices can 

acquire values of 1 or 2. 

Using the Born weak coupling approximation (with 
the hierarchy W > U > V of the model parameters, 
where W is the bandwidth for the upper and lower 
branches of the energy spectrum and © of graphene 
for the case of t 2 =0), we can consider the scattering 
amplitude in the Cooper channel, confining our anal¬ 
ysis to only second-order diagrams in the effective in¬ 
teraction of two electrons with opposite values of the 
momentum and spin and using quantity r(p,k) for it. 
Graphically, this quantity is determined by the sum of 
the diagrams shown in Fig. [TJ Solid lines with light 
(dark) arrows correspond to Greens functions for elec¬ 
trons with opposite values of spin projections (~^)- 
It is well known that the possibility of Gooper pairing 
is determined by the characteristics of the energy spec¬ 
trum near the Fermi level and by the effective interac¬ 
tion of electrons located near the Fermi surface [32] • As¬ 
suming that the chemical potential in doped graphene 
falls into the upper energy band Uik and analyzing the 
conditions for the occurrence of Kohn- Luttinger super¬ 
conductivity, we can consider the situation in which the 
initial and final momenta also belong to the upper sub¬ 
band. This is reflected in Fig. [T] via indices a.\ (upper 
band) and (lower band). 

The first diagram in Fig. [1] corresponds to the ini¬ 
tial interaction of two electrons in the Gooper channel. 
The next (Kohn-Luttinger) diagrams in Fig. [T] describe 
second-order scattering processes, (5r(p,k), and take 
into account the polarization effects of the filled Fermi 
sphere. Two solid lines without arrows in these dia¬ 
grams indicate summation over both spin projections. 
Wavy lines correspond to the initial interaction. Scat¬ 
tering of electrons with identical spin projections cor¬ 


responds only to the intersite contribution. If electrons 
with different spin projections interact, the scattering 
amplitude is determined by the sum of the Hubbard 
and intersite repulsions. Thus, in the presence of the 
short-range Goulomb interaction alone, the correction 
i5r(p, k) to the effective interaction is determined by the 
last exchange diagram only. If the Coulomb interaction 
of electrons at neighboring lattice sites of graphene is 
taken into account, all diagrams in Fig. [T] contribute to 
the renormalized amplitude. 

After the introduction of the analytical expressions 
for the diagrams, we obtain the following analytic ex¬ 
pression for the effective interaction in Fig. [TJ 


r(p, k) = -^ -b -^Re(Mp_krpr’k) + i5r(p, k), 

where 

<5f(p,k) = ^ ^ r5L.^^.(p,q2| - k,pi) X 

xrju;ii(Pi. -p|q2, k)xi,j(q2, Pi) + 

H {ri;*i(p,Pilqi,k) X 

X [r|^yi(qi> -p|Pi> -k) - rl'i.y(qi, -p| - k, p^) 


(16) 

(17) 


+ 


+ra;ij(qi,-p| - k,Pi) X 


(p> Pi |k, qi ) - r'1^.. (p. Pi Iqi , k) |xij (q^, p^ ). 


Here, we have introduced the following notation for gen¬ 
eralized susceptibilities: 


Xij(k,p) 


/(Aik) /(Ajp) 


(18) 


where 

/(a:) = (exp(^) + l)-i 
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— is the Fermi-Dirac distribution function, and ener¬ 
gies ifik are defined by expressions For the sake of 
compactness, we have introduced the notation for the 
combinations of momenta: 

Qi = Pi + P - k, q2 = Pi - p - k. (19) 

Knowing the renormalized expression for the effec¬ 
tive interaction, we can pass to analysis of the condi¬ 
tions for the emergence of superconductivity in the sys¬ 
tem under investigation. It is well known [49] that the 
emergence of Cooper instability can be established from 
analysis of the homogeneous part of the Bethe-Saltpeter 
equation. In this case, the dependence of the scatter¬ 
ing amplitude r(p,k) on momentum k can be factor¬ 
ized, which gives the integral equation for the supercon¬ 
ducting order parameter A(p). After integrating over 
isoenergetic contours, we can reduce the problem of the 
Cooper instability to the eigenvalue problem [551I3DH51] 

where superconducting order parameter A(q) plays the 
role of the eigenvector, and eigenvalues A satisfy the re¬ 
lation ~ ln(rc/IF). In this case, momenta p and q 
lie on the Fermi surface and w_F(q) is the Fermi velocity. 

To solve Eq. dSOj), we write its kernel as the super¬ 
position of eigenfunctions each of which belongs to one 
of irreducible representations of symmetry group C^v of 
the hexagonal lattice. It is well known that this sym¬ 
metry group has six irreducible representations m- 
four ID and two 2D representations. For each represen¬ 
tation, Eq. (EOl) has a solution with its own effective 
coupling constant . We will henceforth use the follow¬ 
ing notation for the classification of the order param¬ 
eter symmetries: representation Ai corresponds to the 
s—wave symmetry type; Bi and B 2 correspond to the 
/—wave symmetry; Ei, to the p + ip—wave symmetry 
type; and E 2 , to the d + id—wave symmetry type. 

For each irreducible representation v, we will seek 
the solution to Eq. (EOj) in the form 

AM(</))=^AMgM(<),), (21) 

m 

where m is the number of the eigenfunction belonging 
to representation v and (j) is the angle determining the 
direction of momentum p relative to the Px axis. The 
explicit form of the orthonormal functions gm\4>) is de¬ 
fined by the expressions 


^i-» gm^W= rr, - ^=cos6to(/), me[0,oo), 

V(1 + Omo)7r 

^ sin (6 to + 6)^, 

^ 1 -^ = ^sin(6m-f 3)/i, (22) 

B 2 ^ = -)= cos (6m + 3)(/), 

Ei^ g^+*^H<^) = ^(Alsin(2m + l)<^ + 

+B cos (2m -I- !)(/), 

E 2 ^ (Asin(2m-f 2)(^-f 

VTT 

+B cos (2m -I- 2)(/). 

Here, subscripts m for the 2D representations Ei and 
E 2 run through the values for which coefficients (2m -I- 1) 
and (2m -I- 2), respectively, are not multiples of three. 

The basis functions satisfy the orthonormality con¬ 
ditions 

27r 

J d(t)g^rl^H(j))gi^\(t)) = S^lsSmn- (23) 

0 

Substituting expression (EU into Eq. (EHI), integrat¬ 
ing with respect to angles, and using condition (l23ll . we 
obtain 

^AMaM=A.AM, (24) 

n 

where 

Ett Ett 

= |q) X 

0 0 

X (25) 

Since Tc ^ IFexp(l/A), each negative eigenvalue 
corresponds to the superconducting phase with the or¬ 
der parameter symmetry v. Generally speaking, the 
expansion of the order parameter (</) in eigenfunc¬ 
tions includes a large number of harmonics; however, the 
main contribution is determined by only some of these 
harmonics. The highest value of the superconducting 
transition temperature corresponds to the modulus of 
the largest value of X^. 

Figure shows the calculated dependencies of the 
effective coupling constant A on carrier concentration n 
for various symmetry types of the superconducting or¬ 
der parameter for the set of parameters t 2 = 0, U = 
2|ti|, and V = 0. It can be seen that for low elec¬ 
tron densities 1 < n < 1.12, in the vicinity of the 
van Hove singularity, the competition occurs between 
the superconducting phases with the /—wave symmetry 
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Fig. 2. Dependences of A on carrier concentration n 
in the graphene monolayer: (a) t 2 = 0,U = 2|ti|, and 
V = 0] (b) t 2 = 0,U = 2|ti|, and V = 0.5|ti|. 


type, whose contribution is determined by the harmon- 

( f \ 1 

ics gm {(f)) = sin (6 to + 3)(j), and the d + id—wave 

symmetry type corresponding to 2D representation E 2 ■ 
In the electron concentration range 1 < n < 1.12, the 
d +id—wave pairing prevails, while for 1.12 < n < n-uH, 
superconductivity with the /—wave symmetry type of 
the order parameter is stabilized. 

It should be noted that to avoid the summation of 
parquet diagrams [SilSI], we do not analyze here the 
electron concentration ranges which are too close to the 
van Hove singularity (Fig. [S]). 

The account of the intersite Coulomb interaction 
considerably affects the competition between supercon¬ 
ducting phases. This can be seen from Fig. which 
shows the A(n) dependences for parameters t 2 = 0,U = 
2|ti|, and V = 0.5|ti|. Comparison with Fig. shows 
that the switching of the intersite Coulomb interaction 
suppresses Cooper pairing in the d + id-wave channel 


did 



n 


Fig. 3. Modification of the electron density of states for 
the graphene monolayer upon switching of the hoppings 
to the next-to-nearest atoms for t 2 = 0 solid curve), 
t 2 = — 0 . 2 / 1 1 (dashed curve), and t 2 = 0.2/i| (dotted 
curve). 


for low electron densities; however, it leads to super¬ 
conductivity with this type of symmetry in the vicinity 
of the van Hove singularity. As a result, the /-wave 
pairing takes place in the electron concentration range 
1 < n < 1.22. It should be noted that this result is in 
qualitative agreement with the results obtained in m- 

The switching of electron hoppings t 2 to the next-to- 
nearest carbon atoms in the graphene monolayer does 
not qualitatively affect the competition between the 
superconducting phases of different symmetry types, 
which is illustrated in Fig. [SB]. Such a behavior of 
the system can be explained by the fact that an account 
of hoppings f 2 > 0 or ^2 < 0 does not noticeably modify 
the density of electron states of the monolayer in the 
range of carrier concentrations between the Dirac point 
and both van Hove singularities (Fig. [3]). However, the 
inclusion of hoppings t 2 leads to an increase in the ab¬ 
solute values of the effective interaction and, hence, to 
higher superconducting transition temperatures in the 
idealized graphene monolayer [58) . 

3. IDEALIZED GRAPHENE BILAYER 

Let us consider an idealized graphene bilayer, as¬ 
suming that two monolayers are arranged in accordance 
with the AB type (i.e., one monolayer is turned through 
60° relative to the other monolayer) (SUED]. We choose 
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Fig. 4 . Crystal structure of the graphene bilayer. Car¬ 
bon atoms Ai and Bi in the lower layer are shown by 
red and black balls; atoms A2 and B2 in the upper layer 
are shown by black and green balls. Intralayer electron 
hoppings are marked by ti and t2', 71, 73, and 74 show 
the interplanar hoppings. 


Here, we have used notation analogous to that for 
Hamiltonian (HD for a monolayer in Section 2. Index 
f = 1,2 in Hamiltonian (I^Sl) denotes the number of the 
monolayer. We assume that one-site energies are iden¬ 
tical (eai = £A 2 = £bi = £b 2 = e)- Interlayer elec¬ 
tron hopping parameters are denoted by 71, 73, 74 (see 
Fig. 2]). The last three terms in Hamiltonian (051) take 
into account the interlayer Coulomb interaction of elec¬ 
trons in atoms Ai and H 2 , Bi and B 2 , and Ai and B 2 ', 
the intensities of these interactions are denoted by Gi, 
G 3 , and G 4 , respectively. 

Passing to the momentum space, it is convenient to 
write Hamiltonian Hq in matrix form: 




E 

k(T 


/ a! 


Ikcr 




^ 2 ka 

^Iko- 

V ^ 2 ko- ) 


Ek 

7i 

tlUk 

V 74Mk 


7i 

Ek 

74 Uk 
tiK 


GUk 74Wk ^ 



74Uk tiUh 


^2kcr 

Ek 73 Mk 


^lk(7 

73 Mk Ek / 


^ ^2kcr 


the arrangement of the sublattices in the layers in such 
a way that the sites from different layers located one 
above another belong to sublattice H, while the remain¬ 
ing sites belong to sublattice B (Fig. 0]). In this case, 
the Hamiltonian of the graphene bilayer in the Wannier 
representation has the form 


H = Ho + Hint, 

^o = (e-m) (e^e+e 


(26) 




,i/CT 


iga 


tl ^ 'X ^\fa-^IJ+S,cr + 02 /(T^ 2 ,/-( 5 ,o- + h.C.) 
fSa- 

— t 2 I alj:^aimcr + + h.C. | 

i-T \«/m» iign)) J 

- 71 + h.C.) - 73 ^{b\g^b 2 ,g+ 5 ,c^ + h.C.) 

fa gSa 

~ 74 ^ ] {^\f^b2J-S,a + 02/(T^1./+'5,'T h.C.), (27) 


fSa 


Hint - ( E 'bft'Hfi w ^ 'Hgt'Hgi 

if ig 


hff-thff I + 


V 


y^,(^ffa'b'f,f+S.a + f'‘2fa'b'2J-S,a 


fSa 


where Sk = fo/k — e, and quantity fit: is defined by ex¬ 
pression ©■ 

Hamiltonian Hq can be diagonalized using the Bo- 
goliubov transformation 

aikcr = Wii (k)aiko^ + Wi2 (k)a2ko- + (30) 

-I- Wi3(k)biita- + Wi4i'k)b2ita- * = 1,2, 3,4, 

It acquires the form 

4 

Hq — Gikalko-Q^ikcr- (31) 

i—1 ktr 

According to the results of ismsu, the interlayer hop¬ 
pings 74 are relatively weak, so it allows us to assume 
that 74 = 0 for convenience of diagonalization of the 
Hamiltonian. In this case, the four-band energy spec¬ 
trum of the graphene bilayer is described by the expres¬ 
sions 

Eik = E ± Y^Ak ± \/Bk — fo/k, (32) 

Ak = ^ ^2a^ -|- 4|6kP + 2|dkP^ , 

Hk = — ^|dk|^(Mk|^ ~ 2a^ + 4|6k|^) + + 4a^|6k|^ + 

+Aab^dl. + 4a6^dk), 

a = 7i, bk = tiUk, dk = 73Wk, 


-b Gi E f^lfan2fa + ^3 XI 

fa gSa 

+ ^4 X {f‘‘lfa'^2,f-S,a + '^2fa'b'lJ+5,a^ ■ 
fSa 


where quantity Uk is defined by expression (|S]). 

Analysis of the conditions for the occurrence of 
(28) Kohn-Luttinger superconductivity in the graphene bi¬ 
layer is carried out in accordance with the general 
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Fig. 5. (a) Energy structure of the graphene bilayer 
near Dirac points and (b) formation of the multisheet 
Fermi contour at weak doping. 



Fig. 6. Dependence of the electron density of states 
for the graphene bilayer per unit cell of one layer on 
the electron concentration for the set of parameters 
t2 = 0, 7i = 0.12|ti|, 73 = 0.1|ti|. 


scheme described in Section 2. We consider the situ¬ 
ation in which, as a result of doping, the chemical po¬ 
tential of the bilayer is in the two upper energy bands 
Elk and E 2 k as shown in Fig.[S^. In this case, the initial 
and final momenta of electrons in the Cooper channel 
also belong to the upper two bands; for this reason, 
indices i and j in the Kohn-Luttinger diagrams for a 
bilayer (see Fig. [1]) acquire the values 1 or 2. Writing 
the analytical expressions for the diagrams, we obtain 
the analytic expression for the effective interaction of 
electrons in the Cooper channel of the graphene bilayer 
in Fig. [U which can subsequently be used for analyzing 
the homogeneous part of the Bethe-Saltpeter equation. 
When solving eigenvalue problem (uni), integration is 
carried out with the allowance for the multisheet na¬ 
ture of isoenergetic contours (Fig. O)). 



Fig. 7. Dependence of A on carrier concentration 
n in the graphene bilayer for the set of parameters 
= 0, t/ = 2, 71 = 0.12, 73 = 0.1, V = 0.5, Gi = 
Gs = G 4 — 0.5 (all parameters are given in the units of 

N). 


Let us now consider the dependences of effective cou¬ 
pling constant A on carrier concentration n for various 
types of symmetry of the superconducting order param¬ 
eter in the graphene bilayer. It should be noted that 
in numerical calculations for the graphene bilayer for 
7 i = 73 = 74 = 0 and Gi = G 3 = G 4 = 0 , we get a lim¬ 
iting transition to the results obtained in Section 2 for 
a graphene monolayer. Figure [7] shows the \{n) depen¬ 
dences determined for the bilayer with the set of param¬ 
eters t 2 = 0,U = 2\ti\, 71 = 0 . 12 |ti|, 73 = 0 . 1 |ti|, 74 = 
0, and F = Gi = G 3 = G 4 = 0.5|ti|. The values of the 
intralayer and interlayer hopping integrals used here are 
close to the values determined in iniiisa for graphite. 
The electron density of states for the graphene bilayer 
for the given set of parameters is shown in Fig. [B] To 
demonstrate the effect of the interlayer Coulomb inter¬ 
action, we chose the maximal possible values of intensity 
Gi, G 3 , and G 4 for which the weak-coupling approxima¬ 
tion is still applicable. Comparison with Fig. [5 ]d shows 
that the allowance for the interlayer interactions does 
not change the domains of superconductivity with the 
f— and d -I- id-wave symmetry types; however, it leads 
to a significant increase in the absolute values of the ef¬ 
fective coupling constant and, hence, to an increase in 
the superconducting transition temperature. 


4. CONCLUSIONS 

We have analyzed the conditions for the emergence 
of superconducting Kohn-Luttinger pairing in systems 
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with a linear dispersion relation using as an example an 
idealized graphene monolayer and bilayer, disregarding 
the Van der Waals potential of the substrate and impu¬ 
rities. The electronic structure of graphene is described 
using the tight binding method in the Shubin-Vonsovsky 
model taking into account not only the Coulomb repul¬ 
sion of electrons on the same carbon atom, but also 
the intersite Coulomb interaction. It is shown that 
the inclusion of the Kohn-Luttinger renormalizations 
up to the second order of perturbation theory inclu¬ 
sively and the allowance for the intersite Coulomb in¬ 
teraction determine to a considerable extent the com¬ 
petition between the superconducting phases with the 
/— and d + id-wave types of the order parameter sym¬ 
metry. They also lead to an increase in the absolute 
values of the effective interaction and, hence, to higher 
superconducting transition temperatures for the ideal¬ 
ized graphene monolayer. 

The results obtained for the graphene monolayer are 
generalized to the case of an idealized graphene bilayer 
consisting of two monolayers interacting via Coulomb 
repulsion between the layers. It is shown that the anal¬ 
ysis of the idealized bilayer system of graphene leads to 
a considerably higher value of the superconducting tran¬ 
sition temperature in the context of the Kohn-Luttinger 
mechanism. 
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